Non-local network-excitation in the Ising antiferromagnet on the anisotropic 
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We study the Ising antiferromagnet on the triangular lattice which are interacting along three 
different directions via two J-bonds and one anisotropic J'-bond with J' > J. Although its finite 
temperature state has long been considered as simply disordered, we find a systematic generation 
of a ferromagnetic J'-bonds named "good defects" each carrying the energy 25j = 2{J' — J). They 
exhibit an eminent correlation which is regarded as a network or a soft lattice. The specific heat 
shows a universal peak at T* ~ Sj and the corresponding energy follows the activation type of 
temperature-dependence with an excitation gap 25 j. Energetics and the crossover between this 
particular two-dimensional state and the other disordered or one- dimensional states are discussed. 
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The Ising model continues to be the starting point of 
understanding the basic concept of magnetism from the 
date of its proposal^]. In the midth of the last cen- 
tury, the exact solutions by OnsagerQ, WannierQ, and 
others 0, successively clarified its nature on the two 
dimensional (2D) lattices. Particularly, the phase transi- 
tion at finite temperature found on the square lattice Q 
remains a prototype description of the conventional phase 
transition and critical phenomena. Meanwhile, still at 
present, there are representative models whose low en- 
ergy properties are studied extensively, some of them are 
related to a series called the frustrated systems [6]. 

Triangular lattice systems are the representative geo- 
metrically frustrated systems which are studied in the 
context of the SU(2) J,,^, anisotropic l3, spin 
models, and in electronic systems as well which yields 
a spin Hquid Mott insulator [isl 14 1. While the ground 



state of these models are quite often disputing, the low 
energy excitations of such frustrated systems are far more 
difficult, where we expect the dense spectrum due to 
the many-body effect. One of the clue to understand 
them is to start from the Ising limit and include quan- 
tum or thermal fluctuations la, 1^1 ■ The present Letter 



actually starts from the exact solutions of the triangular 
Ising antiferromagnet Q which is classically disordered 
throughout the whole temperature range. Revisiting this 
"already well understood" system, we clarify a new as- 
pect of thermal fluctuation in the anisotropic but still 
heavily frustrated region. The Hamiltonian is given as, 



(1) 



where = ±1 and interactions are confined to nearest 
neighbor sites only. Throughout this Letter we focus on 
the case where the interactions along bonds in one di- 
rection Jij = J', which is taken as a unit of energy, are 



larger than in the other two directions, Jij — J, namely 
< J < J' . Hereafter we regard the direction along the 
strong J'-bonds as the chain direction. The ground state 
at J < J' is a combination of one-dimensional (ID) Neel 
order along the strong J' chains and a disorder in the 
remaining two directions of J-bonds where the spin-spin 
correlation along decay in power law. We find that at 
finite temperature, the partial disorder embedded in this 
system couples with the thermal fluctuation and yields 
an exotic state where the ferromagnetic J'-bonds form 
soft non-local quasi-2D network. 

We flrst analyze thermodynamical properties obtained 
by the exact solution^. Figures [T^a) and [T](b) present 
the exact energy e = E/N and the specific heat c = 
de/dT for various values of J ^ 1. There is no disconti- 
nuity in these bulk properties. We find, however, a peak 
or its precursor of the specific heat at low temperatures 
when J > 0.5 J'. The inset of Fig. [TJa) shows an Arrhe- 
nius plot of the dimensionless energy (with respect to the 
ground state energy Cg — — J') in unit of Sj = J' — J as 
a function of dimensionless temperature. We find that 
these curves asymptotically approach an activation-type 
behavior with the energy gap A = 26 j. This character- 
istic energy scale is also found in the universal behavior 
of the specific heat in the inset of Fig. [ijb) as a func- 
tion of normalized temperature, T/Sj, where the peak at 
T/6,j ~ 1 is scaled for all cases of 0.5 < J < 1. 

Let us examine the possible configuration of spins at 
low energy in order to understand the origin of the char- 
acteristic energy scale A. We start from the disordered 
ground state as shown in Fig. [2ja). Here, each chain 
has a Neel order. Inter-chain interactions J are cancel 
out on the whole irrespective of a relative phase of an- 
tiferromagnetic spin configurations between the chains, 
as shown in the right panels of Fig. [2fa), and thus the 
ground state has degeneracy. An excited state is gener- 
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FIG. 1: Panels (a) and (b) are the energy per site, e — E/N, 
and the specific heat, c = de/dT, as a function of temperature 
T for various choices of J ^ J' = 1. The inset of panel (a) 
gives an Arrhenius plot for (e — eg)/Sj as a function of scaled 
temperature T/Sj. Solid line gives the activation-function 
with A = 25 J. Inset of panel (b) gives the replot of c as a 
function of ksT/Sj. Panel (c) shows the direct comparison 
of energy the Monte Carlo simulation (marks) and the exact 
solution (solid lines) at J = 1 and 0.95. 



ated by introducing a ferromagnetic configuration on a 
J'-bond as illustrated in Fig. [DJb), which we hereafter 
refer to as a "defect". Each defect takes an excitation 
energy of 2J' along the chain direction, while between 
the chains it creates an uncanceled inter-chain interac- 
tion energy. In order to examine the energetics of the 
defects in detail, we decompose the lattice into plaque- 
ttes depicted in Fig. (He), and count the energy of them 
by summing up the energy of J'- and J-bonds within the 
plaquette. Here, we take half of the energy of the four J- 
bonds on the edges of the plaquettes which are shared by 
the neighbors. Plaquettes are energetically classified into 
four types. Those without defects are assigned (I) with 
an energy — J' , where the interchain interactions are in- 
deed canceled out. Each plaquette categorized in (II) has 
a defect which gains the interchain bond energy, and the 
energy amounts to J' — 2 J. Those classified as (III) and 
(IV) have higher total energies, J' and J' + 2 J, respec- 
tively. Thus, in the region J < J' , plaquettes (I) have the 
lowest energy, while plaquettes (II) have the first excited 
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FIG. 2: Configurations of (a) ground state and (b) a single 
defect state. Black and white circles represent the up and 
down spins, respectively. The J'-bond which connects iden- 
tical spins are called "defects", (c) Classification of the unit 
plaquettes according to their energy. See main text for details. 



state energy which we hereafter call as "good" defects, 
and the plaquettes (III) and (IV) will be contrastingly 
referred to as "bad" defects. The excitation energy of 
defect (II) actually correspond to the characteristic en- 
ergy scale A = 26 j in the previous exact solution. If Sj is 
sufficiently small, the low energy excitations by the bad 
defects can be negligible and the state with only "good" 
defects may dominate. 

In order to examine to what extent the good defects 
dominate at low temperatures, we perform the classical 
Monte Carlo simulation (MC). Calculations are carried 
out mainly on a iV = iV^; x Ny = 50 x 100 site cluster up 
to the observation time of 5000 Monte Carlo steps (MCS) 
after the relaxation (which takes less than 500 MCS). We 
take about 200 independent runs for averaging. For each 
snap shot we separately count the number of good defects 
(Ng) from the total number of defects (Nd). We confirm 
that the energy of MC almost exactly coincides with that 
of the exact solutions as in the inset of Fig. [ijc) by the 
dotted (MC) and solid(exact) lines. 

We first present the T-dependence of = Nd/N and 
Ug = Ng/N at J = 0.95 in Fig. [3{a) for various clus- 
ter sizes. We confirm that size effects are negligible 
here. Both Ud and Ug decreases exponentially as T is 
lowered. At higher temperature, the number of bad de- 
fects, rifc = nd — Ug, increases rapidly. A distinct fea- 
ture is found in the T-dependence of Ug/ud in Fig. E^b) 
at around the particular degree of anisotropy, J = 0.5; 
when J > 0.5, the ratio Ug/ud increases monotonically 
toward unity with decreasing temperature, whereas at 
J < 0.5, a peak exists at T = T,n- The behavior of 
Ug/ud is better understood in its contour plot on the J- 
T diagram given in Fig. [3l^c). Now, we plot Tm in a J-T 
diagram in Fig. [S^d). In the shaded region the good de- 
fects are energetically chosen to be thermally activated. 
Here, fc^T = 2J which is the energy difference between 
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FIG. 3: The results of the Monte Carlo simulations at J' = 
1. (a) Density of defect (II- IV), rid, and good defects, Tig, 
for different sizes of clusters as a function of temperature at 
J' = 1. (b) Ug/ud as a function of temperature, (c) Contour 
plot of ng/ud by J and fcsT. (d) J-T diagram which depict 
the dimensional crossover. 



bad (III) and good (II) defects, approximately gives the 
crossover line which characterizes the temperature that 
the bad defect is activated. Above this crossover line, 
the difference between the good and bad defects no more 
makes sense. Rather, any defects are activated randomly 
as the temperature is increased. 

Let us next discuss the nature of the good defect states 
in detail. We find in the MC results that defect-defect 
correlation decays quickly toward zero after a few chain 
distance in the interchain-direction (not shown). Then, 
the defects seemingly distribute randomly in real space, 
which is however, not the case. Instead, an distinguished 
correlation arises in the positions of good defects as a 
consequence of eliminating bad defects. 

One finds the following configuration rule of the good 
defect state; to place the good defects exclusively, those 
on two neighboring chains must always be located alter- 
natively along the chain. We indeed confirm in the MC 
that all the snap shots with Ug/nd = 1 strictly follows 
this rule. The representative good defect configuration 
is given in Fig.|4l^a). Here, all the defects are good ones, 
and if we neglect the background spins as in Fig. Hfb), 
the good defects form soft lattice structure. Here, we 
mean by "soft" that the absolute location of defects do 
not matter, i.e., defects are shifted in the chain direction 
without breaking the above rule only if it does not go 
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FIG. 4; (a) Representative configuration of the "good defect 
state" with Ugjud = 1. Arrows indicate the two sites identi- 
cal under the periodic boundary condition along the chains. 

(b) Schematic "soft lattice structure" of good defects corre- 
sponding to (a). Panels (c), (d) and (e) describes the multiple 
defect configuration embedded in the antiferromagnetic back- 
ground. See main text for details. The lower panels of each 

(c) -(e) show the equivalent case when the defects are shifted 
and squeezed along the chain, (f) The cases where the two 
defects on the center chain are depleted. 



over the defects on the neighboring chains. 

We now explain the energetics behind the above rule 
step by step. It is noteworthy that the relative phase of 
the Neel state between the neighboring chains essentially 
determines the energy of a defect when it is embedded 
as shown in Fig. [D^b). We start from this figure noting 
that the position of the defect can be shifted translation- 
ally along the chain without changing its energy. This 
is because such a shift by unit length flips all the spins 
in a defect-plaquette, in which case a good (bad) defect 
remains good(bad) with the energy unchanged. Next, we 
consider multiple defects on a chain. Note that if we re- 
strict ourselves to the periodic boundary condition with 
even number of sites along the chain, the defects should 
always appear as pairs. Figure |D^c) shows two defects on 
a single chain. A pair of defects on isolated chain is al- 
ways a combination of good and bad ones. This is easily 
understood if the defects are shifted and placed next to 
each other as in the lower panel of Fig.|4l^c). Thus, as long 
as we consider pairs of defects on an isolated chain the 
bad defects cannot be excluded. Now, in order to make 
both defects good, one need to twist the phase of antifer- 
romagnetic configuration on the neighboring chain. This 
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corresponds to inserting one defect into the neighboring 
chain between the previous two as in Fig. [4jd), and the 
introduced defect also turns out to be a good one. If we 
insert one more defect as in Fig. Sfe), we again cannot 
exclude bad defects, since two consecutive defects bring 
the phase of the antiferromagnetic configuration back. 
The picture is easily extended to multiple chains, and 
the configuration pattern in Fig.lDJd) is extended to that 
of Fig. HJa). This straightforwardly leads to the "good 
defect rule" we mentioned earlier. 

The remaining issue is the crossover at J = 0.5, which 
is related to the stability of the soft lattice state. Con- 
sider depicting a pair of defects on one of the chains in 
this state. The equivalent case where a good defects are 
packed closest together is given in Fig. |DJf). The deple- 
tion corresponds to replacing two (Il)-plaquettes on the 
chain by (I) and two (Il)-plaquettes in the neighboring 
chains in between the depleted two by (III), and the net 
energy gain becomes dE = 2 x {2J — J'). At J < 0.5J', 
we have dE < and the soft lattice is energetically un- 
stable. Then, the isolated defects shown in Fig. Ufc) can 
be realized at low temperatures, namely at T < and 
J < 0.5 in Fig. [3jd). Note that this configuration of 
defects corresponds to the classical version of fractional 
charges in the systems with quantum fluctuation [l6i|. 

In the presence of the soft lattice structure, density 
fluctuations of defects show characteristic behaviors; it is 
disfavored to change numbers of defects independently on 
each chain since such fluctuations destroy the soft lattice 
structure. Instead, defects must be inserted/depleted si- 
multaneously on all chains in such a way that the good 
defect rule is maintained. Actually, in MC when an ad- 
ditional defect is introduced on a certain chain, it con- 
tinuously spreads to the neighboring chain to form a soft 
domain wall or a soft dislocation line in the interchain 
direction to make all chains finally equivalent. This na- 
ture of the defect density fluctuations changes drastically 
once we go over the crossover line atJ = 0.5. AtJ<0.5, 
it is energetically allowed to change the defect densities 
on each individual chain independently. Namely, we may 
regard that the defect density fluctuation have ID char- 
acters, similar to that found in the limit J — 0. 

The J-T diagram in Fig. ^d) is thus well classified 
by the dimensional crossover lines. In the whole region 
of ^ J < 1, spins in the ground state exhibit a ID 
Neel configurations along the chain, and have degener- 
acy in such a way that all the spins on a single chain 
can be flipped simultaneously without any cost of energy. 
Namely, the interchain coupling J is irrelevant. However, 
once we enter finite temperature region at J > 0.5 where 
the 2D soft lattice network of good defects is present, 
such ID fiuctuation disappears. This is because if we 
flip all the spins simultaneously on one chain in the soft 
lattice state, the good defects on this chain and on the 
neighboring two chains are all transformed to bad ones, 
and the energy increases. In other words, the relative 



phases between the chains are locked in forming the soft 
lattice state. Then, instead of the above mentioned ID 
fluctuations, the 2D defect density fluctuations appear, 
which allows for the change of defect number on all chains 
simultaneously. The ground state at J = 0.5J' is thus 
a singular point where a dimensional crossover between 
T = and T > emerges. 

To summarize we studied the triangular Ising antifer- 
romagnet interacting with one strong anisotropic J'- and 
two weaker J-bonds. The system systematically gener- 
ates ferromagnetic J'-bonds (defects) under thermal fluc- 
tuation, which exhibit an eminant non-local correlation 
throughout the system by conserving their number and 
by keeping their relative locations. This state which is 
regarded as two-dimensional soft network structure, can- 
not be described by the conventional " long range order" , 
since it is neither characterized by the spin-spin nor the 
defect-defect correlation functions. The present flndings 
will provide clue to understand the low-energy excita- 
tion of the anisotropic triangular antiferromagnet under 
fluctuations from the Ising limit |15l. M, Il3- As far as 
we know, such novel phenomena has never been found, 
and shall provide a new example in exploring the scheme 
beyond the conventional statistical or condensed matter 
physics. 
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